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Abstract. The dynamics of steps on crystal surfaces is considered. In general, the 
meandering of the steps obeys a subdiffusive behaviour. The characteristic asymptotic 
time laws depend on the microscopic mechanism for detachment and attachment of 
the atoms at the steps. The three limiting cases of step-edge diffusion, evaporation- 
condensation and terrace diffusion are studied in the framework of Langevin descrip- 
tions and by Monte Carlo simulations. 



1 Introduction 

Surface dynamics has attracted much interest in recent years, partly due to 
progress in experimental techniques to measure atomic processes at surfaces. 
A variety of phenomena has been studied, both experimentally and theoreti- 
cally, including fluctuations of isolated and bunches of steps, island motion, and 
thermal relaxation of sinusoidal gratings (Lancon and Villain (1990), Pimpinelli 
et al. (1993), Williams (1994), Selke and Duxbury (1995), Hager and Spohn 
(1995), Khare et al. (1996), Bonzel and Surnev (1997), Weeks et al. (1997), Khare 
and Einstein (1998)). 

In this contribution, we shall draw attention to the (sub-) diffusive dynam- 
ics of the meandering of steps of monoatomic height. We shall emphasize the 
role of different atomic mechanisms leading to distinct limiting cases for the 
asymptotic time laws. In particular, (i) evaporation-condensation (EC), i.e. un- 
correlated attachment and detachment of atoms at the step edge, (it) periphery 
diffusion (PD) along the step edge, and (iii) terrace diffusion (TD) will be dis- 
criminated (these three cases correspond to evaporation-condensation, surface 
diffusion and volume diffusion in the pioneering work of Mullins (1959) on ther- 
mal flattening of crystal surfaces). Results of Langevin descriptions and of Monte 
Carlo simulations will be presented. 



2 Step Fluctuations 



Let us consider a surface with an isolated step of monoatomic height. The step 
position, at time t and step site x, is specified by h{x, t) (in units of the lattice 
spacing, i.e., h = ... — 1,0,1,...; likewise x = 1,2,..., L, where L is the length 
of the step). For instance, h{x,t) — describes a perfectly straight step. Step 
fluctuations may be quantified by the correlation function 

G{t) {h{x, t) ~ h{x, 0)f >x.ensemble (1) 
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where the brackets denote spatial (along the step edge, x) and ensemble averages. 
The initial step position, h{x, 0), may be chosen in various ways, for instance, it 
may be straight or thermally equilibrated. 

At fixed time, the step stiffness may be characterised by the spatial correla- 
tion function 

F{Ax) {h{xo, t) - h{xi,t)f >ensemble (2) 

where Ax is the distance between xq and xi. Choosing xq and a^i in the center 
and at the end of the step, respectively, F then measures the width of the step 
(and will be denoted by w). Usually, the width of the step is limited by and 
related to its length, L. 

If each step position would be completely free to wander back and forth with 
the same probability at any time (i.e. each step position would be a random 
walker), then G(t) would obviously obey a diffusive behaviour, G{t) oc t. This 
situation may be realized, e.g.. at infinite temperature. However, in general, the 
motion of each step position is hindered due to the other parts of the step. 
Indeed, the step meandering is typically subdiffusivc, with a simple power law 
at late times, G{t) oc i", where the exponent a < 1. The value of a depends 
on the atomic process driving the step fluctuations, and a may be determined 
theoretically or in experiments. 

In the theoretical approaches, assumptions are made on the transition rates 
and correlations for the attachment and detachment of the atoms at the step 
edges as well as on the boundary conditions at the ends of the step (for exam- 
ple, periodic or pinned). In particular, Langevin descriptions and Monte Carlo 
simulations have been applied. 

The basic equation of the Langevin theory for step dynamics has the form 
(Khare and Einstein (1998), Blagojevic and Duxbury (1997)) 

dh{x,t)/dt = J{h{x,t))+r){x,t) (3) 

where the functional J describes the reduction of the step free energy, and r] is 
the noise term. Both quantities depend on the microscopic processes driving the 
step fluctuations. Usually, h and x are assumed to be continuous variables. The 
Monte Carlo simulations (for an introduction, see Binder (1992)) are performed 
for microscopic models, specifying the interactions between the surface (or step) 
atoms and the transition rates for the motion of the atoms at the surface, as 
shown in the following. 

2.1 Evaporation-Condensation 

Let us asstune that the step fluctuations are due to uncorrelated attachment 
and detachment of atoms at the step edges (as may be realized when the surface 
dynamics is governed by exchange of atoms with the surrounding vapour: kinetics 
of evaporation and condensation, EC). 
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In a Monte Carlo simulation, the step may now be mimicked by a one- 
dimensional SOS model 



where the sum runs over all nearest neighbour pairs (x, x') of step sites; J is the 
kink energy. A resulting step configuration is depicted in Fig. 1. 



Fig. 1. Typical step configuration based on the one-dimensional SOS model 
{ksT / J — 1.2), using EC. Dark atoms constitute the upper terrace, light ones the 
lower terrace, with the step sites along the horizontal direction. Only part of the step 
is shown 

The EC kinetics corresponds, in the simulations, to Glauber dynamics. The 
transition rate for attachment or detachment of an atom at a randomly chosen 
step site is given by the Boltzmann factor g-^^/kBT^ with 5E being the energy 
change needed for the possible attachment or detachment, fcs is the Boltzmann 
constant, and T the temperature. Choosing the initial step positions h{x,t = 0), 
G{t) is computed by averaging over a large number, N, of samples. The results 
depend on the temperature, the number of active step sites, L, the boundary 
conditions at the ends of the step, and the initial step configuration. In particular, 
we monitored G{t) and the step width w for steps of length L ranging from 4 
to 200 at various temperatures {kBT/J going from 0.5 to 4.0), pinning the 
step positions at the ends, h{l,t) = hi and h{L + 2,t) — h2 (usually, we took 
/ii = /i2 = 0), or using periodic boundary conditions (pbc) at the ends of the 
steps. The time t is measured in Monte Carlo steps per site (MCS), i.e. one time 
unit has elapsed after L attempts to change the step position by one. Averages 
have been taken over typically N = 10^ to 10^ samples. 

At very short times, the step fluctuations are nearly diffusive, with G{t) oc t, 
due to the random events at (almost) equivalent step sites. However, caused 
by the rigidity of the step, deviations show up soon. They may be monitored 




(4) 
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conveniently by the effective exponent 

aeff{t) = dln{G{t))/dln{t) (5) 
If G{t) would obey a simple power law, the exponent would be a constant, 



1.2 




time (MCS) 



Fig. 2. Effective exponent of tlie step fluctuations, ag//, versus time, in MCS, for EC 
at ksT/J = 1.2, using periodic boundary conditions for initially straight steps (step 
length L= 10, 40, and 200), or with thormalizod steps pinned at the ends, h\ = h2 = 0, 
Z/ = 10 and 200. The solid line refers to Oe// = 1/2 



As depicted in Fig. 2, a^ff tends to decrease with time (apart from the be- 
haviour of initially straight steps, where the exponent first increases slightly to 
values larger than one before getting smaller). The rate of the decrease in a^ff 
depends non-monotonically on temperature, being fastest at some intermediate 
values, say, fesT « J. It is higher when using thermally equilibrated initial step 
positions h{x,t = 0), compared to straight steps. For sufficiently long steps the 
effective exponent seems to approach, at large times, a = 1/2, the asymptotic 
value obtained from Langevin descriptions for EC kinetics of indefinitely ex- 
tended steps (Edwards and Wilkinson (1982), Lipowsky (1985), Abraham and 
Upton (1989), Stauffer and Landau (1989), Bartelt et al. (1994), Blagojevic 
and Duxbury (1997)), and observed experimentally, e.g., on a Au(llO) surface 
(Kuipers et al. (1993), van Beijeren et al. (1983)). 

However, as illustrated in Fig. 2, the finite length of the step leads to different 
asymptotics, depending on the boundary conditions. For pbc, fle// may get close 
to 1/2, before it rises at later times when the step width w saturates. The time t^, 
at which the step width has reached its saturation value depends on the length 
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of the step (exact calculations for a simplified model suggest tyj oc (Abraham 
and Upton (1989)), in reasonable agreement with our simulations). Eventually, 
the step behaves effectively like a single random walker with a reduced diffusion 
coefficient, i.e. fle// — ^ 1- 

In the pinned case, the step width also saturates, w oc L^/^ for both boundary 
conditions (Blagojcvic and Duxbury (1997)). Obviously, the motion of the step 
is now hindered because the step is fixed at the ends. Therefore, G{t) remains 
bounded and Oe// tends towards zero. 

2.2 Periphery Diffusion 

We now suppose that the step fluctuations arc caused by direct hops of step 
atoms between neighbouring step sites (in nature, the hops correspond to the 
motion of individual atoms or they may result from a particle exchange mech- 
anism). The step or periphery diffusion (PD) may be simulated again for a 
one dimensional SOS model, equation (4), applying now Kawasaki dynamics 
(Binder (1992)). The transition rates for diffusion of a randomly chosen step 
atom to a neighbouring site are given by the Boltzmann factor of the energy 
change needed for the hop. 

The Monte Carlo simulations for PD are done in complete analogy to those 
for EC, studying the impact of the step length, the initial step configuration, 
temperature, and the boundary conditions. Most runs were done at ksT/J = 
3.0. Averages were taken over N = 10^ to 10^ samples. The step length L ranged 
from 6 to 200. 

At very short times, the step fluctuations are nearly diffusive, before the 
rigidity of the step starts to slow down their growth. The corresponding rate of 
decrease in the effective exponent Ueff of the correlation function G{t) depends 
again in a non-monotonic fashion on temperature, cieff seems to approach the 
value 1/4 for indefinitely long steps, see Fig. 3, in accordance with predictions 
from Langevin descriptions (Khare and Einstein (1998), Blagojevic and Dtixbury 
(1997)) and experimental data for steps on, e.g., Cu surfaces (Giesen-Seibert et 
al. (1993)). 

However, the finite step length L modifies that behaviour at late times, re- 
flecting also the boundary conditions at the ends of the step, as illustrated in 
Fig. 3. Pinning of the steps limits the fluctuations, with aeff decreasing eventu- 
ally to zero when the width of the step has acquired its saturation value. For 
periodic boundary conditions and short steps, the effective exponent shows a 
rather peculiar time dependence, as depicted in Fig. 3. The step fluctuations 
G{t) are bounded, because the average step position is, by definition of PD, 
conserved, in contrast to the EC case where the entire step of finite width may 
move like a random walker at late times. Accordingly, Ueff goes to zero. Note 
the pronounced maximum in Ueff for the short step (L =10, at about 100 MCS, 
see Fig. 3). This feature seems to be typical for longer steps as well (e.g., for 
L = 20, the maximum shows up at about 1600 MCS). The faster growth in the 
step fluctuations might be related to a modulation in the step positions. Then 
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the amplitudes of the 'hill' and 'valley' parts of the fluctuations may increase, 
satisfying the above conservation requirement, until finally these excursions get 
bounded as well. This aspect, which we also noticed for pinned steps, deserves 
further studies. 



1.2 
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Fig. 3. Effective exponent of the step fluctuations, o-e//, versus time, in MCS, for PD 
at ksT / J = 3.0, using periodic boundary conditions for initially straight steps (step 
length i= 10 and 200), or with thermalized steps pinned at the ends, hi = = 0, 
L = 200. The solid line refers to Oe// = 1/4 

In addition, we studied the situation where the step is pinned at some angle, 
i.e., the ends of the steps are at different positions, hi, and /12. The dynam- 
ics seems to be quite similar to that of stcip where the pinning introduces no 
additional kinks {hi = /12), at least at the fairly high temperature which we 
considered, UbT/J = 3. A more detailed analysis seems to be desirable. 

2.3 Terrace Diffusion 

In general, the step dynamics on crystal surfaces is expected to be rather com- 
plex, because atoms detaching from a step may diffuse (anisotropically and in 
external potentials) on the neighbouring lower and upper terraces, colliding with 
other atoms on the terraces. In addition, exchange mechanisms may introduce 
correlations of different ranges between atoms detaching from and attaching at 
the steps. There have been a few attempts to study the complexity of these 
processes in the framework of Langevin descriptions. A rather appealing pic- 
ture has been suggested by Duxbury and Blagojevic (1997), by expressing the 
functional J of the Langevin equation in terms of the chemical potential at the 
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step sites and a 'diffusion kernel' P{1), which denotes the probabihty of an atom 
being originally one lattice spacing in front of the step site x to attach at step 
position X + I, after having undergone an appropriate walk on the terrace. The 
approach has been critically discussed by Einstein and Kharc (1998). In any 
event, it poses the well defined and interesting problem to calculate the diffusion 
kernel for various scenarios. 

For EC, P(l) would be obviously constant, i.e. the wandering atom could 
attach with the same probability at each step site. Indeed, one obtains, via 
equation (3), that the step fluctuations obey, at late times, the power-law G{t) oc 
t^l'^. For PD, P{1) vanishes for I > 1. yielding the asymptotic exponent a = 1/4. 
Perfect terrace diffusion (TD), where the atom executes a standard random walk 
on the terrace in front of a straight step until it attaches at that reference step, 
is described by a diffusion kernel, at large distance /, of the form P{1) oc 
This decay law leads then to subdiffusive step fluctuations with G(t) oc t^^^, 
being intermediate between EC and PD (Blagojevic and Duxbury (1997)). 

The previous calculations of P{1) for TD had been done in the continuum 
limit for the step position h{x,t). We extended those calculations in several ways. 
P{1) was determined by an exact numerical enumeration of the probability to 
visit a given, discrete site of the terrace (Majid et al. (1984)), which is supposed 
to be a square lattice. In doing that, we are dealing with a lattice of finite extent, 
bordered by a step being parallel and opposite (at a distance d) to the reference 
step and two steps starting at the ends and being perpendicular to the reference 
step of length L. The three boundary steps may either reflect or absorb the 
wandering adatom. The boundary conditions allow one to study, for instance, in 
which way the step dynamics changes because of the presence of a neighbouring 
step with a large (reflecting) or small (absorbing) Schwoebel-Ehrlich barrier. 

In particular, we considered an external potential, V, of the form 

V = A/y' (6) 

where y measures the distance of the diffusing adatom from the reference step; 
A (in units of ksT) is the strength of the potential originating from elastic or 
dipolar interactions of the adatom with the reference step (and its opposite step). 
Thence, the hopping probabitities for the wandering adatom are different for 
moves towards and away from the steps. Assuming straight steps, the resulting 
diffusion kernel P{1) is shown to display the same limiting behaviour as in the 
potential-free (A = 0) case, as illustrated in Fig. 4. In particular, for large 
distances d between the reference and the opposite step, P{1) falls off as 1/P 
for I » 1. Only the prefactor of the power-law depends on the strength A of 
the potential, increasing exponentially with A (at least for small values of A). 
Accordingly, the interaction of the adatom with the step may be expected to 
be irrelevant for TD, with G{t) oc t^^^. For smaller distances d, P{1) decreases 
exponentially with I. The exponential decay, due to reflection or adsorption at 
the neighbouring step, leads to a slowing down of the step fluctuations, being 
now in the limiting case of PD, G{t) oc (Blagojevic and Duxbury (1997)). 
Of course, it would only hold if there is a reflecting (infinite) Schwoebel-Ehrlich 
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barrier, because otherwise the neighbouring (opposite) step would be a source 
of adatoms which could attach at the reference step in an uncorrelated manner, 
giving rise to the fast step fluctuations of EC, G{t) oc t^^"^. 
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Fig. 4. Effective exponent of the diffusion kernel, dlnP(l) / dlnl, versus return distance 
I, for a wandering adatom without (.4 = 0) or in an external potential, V = A/y^ . 
The opposite step is reflecting, while the steps at the side are absorbing. The length 
of the reference step is 1200, while the distance to the opposite step is 600. Note that 
the data are, on the scale of the figure, indistiguishable for large return distances. The 
plateau corresponds to the 1/i^-law discussed in the text. The following decay reflects 
the presence of the opposite step 

Furthermore, wc studied the situation, where the reference step is rough. 
The asymptotics of P{1) seems to be unaffected by this modification as well 
(Bisani (1998)), e.g., if the influence of the opposite step is sufficiently strong, 
we recover the same exponential decay in P{1) as for straight steps, with the 
prefactor refiecting the roughness. 

The findings suggest that the corresponding time dependence of the step 
fluctuations at late stages, G{t) oc t^/^, is rather robust against including elastic 
or dipolar interactions of the diffusing adatom with the step and also against 
the roughness of the step. Different time laws, with a deviating from 1/3, may, 
however, result from a competition between terrace and periphery diffusion with 
various sticking coefficients (Khare and Einstein (1998)) or collisions between 
adatoms diffusing on the terrace. The topic may be studied in simulations using, 
for instance, the kinetic Monte Carlo approach with suitable activation energies 
(or using standard Monte Carlo techniques for a two-dimensional SOS model 
with a step (Bartelt et al. (1994)), applying Kawasaki dynamics). Such simula- 
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tions have been done in particular for a model of steps of monoatomic height on 
Ag(llO) (Stebens (1998)), motivated by recent experiments (Li et al. (1996)). 
The terrace diffusion is extremely anisotropic. The steps may become quickly 
very fuzzy, and they arc no longer uniquely defined, so that rather involved 
analyses are needed to identify the characteristic time-law for the step fluctua- 
tions. 

3 Summary 

We considered three limiting cases for step fluctuations at crystal surfaces, where 
the attachment and detachment of atoms at the steps is (i) uncorrelated, EC, (ii) 
due to local moves along the surface edge, PD, or (iii) mediated by random walks 
of adatoms on the neighbouring terraces, TD. In all three cases, the fluctuations 
of an indeflnitely long step are expected to grow subdiffusively at large times, 
G{t) oc t"", with a~ 1/2, 1/4 and 1/3, respectively. 

The flrst two cases have been studied by doing Monte Carlo simulations 
for one-dimensional SOS models. We conflrmed the generic time laws for step 
fluctuations at late stages, with a — 1/2 for EC, and 1/4 for PD. However, due to 
the flnite length of the step and boundary effects at the ends of the step, various 
other typical scenarios are encountered, including the crossover from subdiffusive 
to diffusive behaviour, for EC and periodic boundary conditions. 

The case of terrace diffusion has been discussed in the framework of a recent 
Langevin description. The main quantity is the diffusion kernel, describing the 
return probability of an adatom wandering on the terrace to a step site. Its 
asymptotics, which is supposed to determine the exponent a (being 1/3 for 
perfect terrace diffusion), is found to be robust against realistic interactions of 
the adatom with the step as well as against roughening of the step. 

Situations in which deviations from these limiting cases may be possible have 
been mentioned. 
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